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^The  scattering  cf  HF  radio  waves  by  a  spherical  electron  cloud  in  the  pres¬ 


ence  of  a  magnetic  field  has  been  obtained  in  ray  optics  approximation  for  a 
Gaussian  distribution  of  electron  density.  The  effect  of  a  magnetic  field  is  small 
at  30  MHz,  but  substantial  at  3  MHz.  The  scattering  in  the  forward  direction  is 
independent  of  critical  radius  and  magnetic  field.  The  displacement  of  the  back- 
scattered  ray  from  the  axis  of  symmetry  of  the  cloud  shows  a  strong  dependence 
upon  field  angle  with  maximum  effects  at  ±45°.  The  cross  section  for  the  ordin¬ 
ary  ray  is  strongly  dependent  on  field  angle  with  a  maximum  near  90°,  the  effect 
decreasing  as  the  critical  radius  increases.  The  extraordinary  ray  shows  a 
minimum  cross  section  near  90°,  with  much  less  dependence  on  field  angle.  The 
maximum  scattering  angle  for  underdense  spheres  exhibits  a  moderate  depend¬ 
ence  upon  field  angle  for  the  ordinary  ray,  and  virtually  none  for  the  extraordi¬ 
nary.  It  is  shown  that  the  backscatter  data  for  the  extraordinary  ray  may  be 
analyzed  in  terms  of  an  effective  critical  radius  which  depends  upon  the  gyro- 
maghetic  ratio,  with  a  small  correction  for  the  effect  of  field  angle. 

®The  objective  of  this  report  is  to  provide  knowledge  and  understanding  of  the 
effect  of  the  earth's  magnetic  field  upon  the  propagation  of  electromagnetic  waves 
through  the  ionosphere,  ionospheric  irregularities,  and  artificial  electron  clouds. 
It  is  thus  directed  to  the  furnishing  of  data,  together  with  criteria,  for  the  solu¬ 
tion  of  Air  Force  problems  of  discrimination,  detection,  and  communication. 
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Abstract 

The  scattering  of  HF  radio  waves  by  a  spherical  electron  cloud  in  the  presence 
of  a  magnetic  field  has  been  obtained  in  ray  optics  approximation  for  a  Gaussian 
distribution  of  electron  density.  The  effect  of  a  magnetic  field  is  small  at  30  MHz, 
but  substantial  at  3  MHz.  The  scattering  in  the  forward  direction  is  independent  of 
critical  radius  and  magnetic  field.  The  displacement  of  the  backscattered  ray 
from  the  axis  of  symmetry  of  the  cloud  shows  a  strong  dependence  upon  field  angle 
with  maximum  effects  at  ±45°.  The  cross  section  for  the  ordinary  ray  is  strongly 
dependent  on  field  angle  with  a  maximum  near  90°,  the  effect  decreasing  as  the 
critical  radius  increases.  The  extraordinary  ray  shows  a  minimum  cross  section 
near  90°,  with  much  less  dependence  on  field  angle.  The  maximum  scattering 
angle  for  underdense  spheres  exhibits  a  moderate  dependence  upon  field  angle  for 
the  ordinary  ray,  and  virtually  none  for  the  extraordinary.  Jt  is  shown  that  the 
backscatter  data  for  the  extraordinary  ray  may  be  analyzed  in  terms  of  an  effec¬ 
tive  critical  radius  which  depends  upon  the  gyromagnetic  ratio,  with  a  small  cor¬ 
rection  for  the  effect  of  field  angle. 

The  objective  of  this  report  is  to  provide  knowledge  and  understanding  of  the 
effect  of  the  earth's  magnetic  field  upon  the  propagation  of  electromagnetic  waves 
through  the  ionosphere,  ionospheric  irregularities,  and  artificial  electron  clouds. 

It  is  thus  directed  to  the  furnishing  of  data,  together  with  criteria,  for  the  solution 
of  Air  Force  problems  of  discrimination,  detection,  and  communication. 
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Scattering  of  HF  Radio  Waves  by  a  Spherical  Electron 
Cloud  in  the  Presence  of  a  Magnetic  Field 


1.  INTRODUCTION 

During  the  past  several  years,  the  Air  Force  Cambridge  Research  Laborato¬ 
ries  have  been  engaged  in  studies  of  the  upper  atmosphere  through  the  release  of 
artificial  ionized  clouds.  A  part  of  this  study  involves  the  use  of  radar  to  obtain 
the  electron  density  distribution  and  its  maximum  value  within  the  cloud.  To  aid 
in  the  interpretation  and  correlation  of  data,  a  systematic  study  is  being  conducted 
on  the  scattering  of  high-frequency  radio  waves  by  electron  clouds.  Since  the  cloud 
size  is  large  in  comparison  to  the  wavelength,  the  investigation  is  being  made  with¬ 
in  the  framework  of  geometrical  optics. 

For  the  case  of  zero  magnetic  field,  the  scattering  amplitude  may  be  obtained 
as  a  function  of  scattering  angle  by  evaluation  of  the  scattering  integral  (Rohm, 

1951).  This  procedure  has  been  utilized  in  order  to  obtain  the  scattering  from  a 
spherical  electron  cloud  with  a  Gaussian  distribution  of  electron  density  in  the  radi¬ 
al  direction  (Klein  and  Mabee,  1969).  When  the  effect  of  the  magnetic  field  is  taken 
into  account,  the  index  of  refraction  becomes  a  function  of  ray  direction  as  well  as 
of  position.  As  a  consequence,  die  scattering  amplitude  cannot  be  obtained  directly. 
Instead,  the  ray  paths  must  first  be  computed  by  a  suitable  technique  and  the  scat¬ 
tering  then  obtained  from  the  angular  displacement  of  neighboring  rays  as  a  func¬ 
tion  of  impact  parameter. 
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To  obtain  the  ray  paths,  we  have  utilized  the  method  developed  by  Haselgrove 
(1955),  which  expresses  the  ray  propagation  in  terms  of  a  set  of  first-order  differ¬ 
ential  equations  in  space  and  time.  The  form  of  these  equations  is  particularly 
amenable  to  solution  by  high-speed  digital  computers.  As  an  initial  phase  of  the  in¬ 
vestigation,  it  is  assumed  that  the  cloud  is  spherical,  with  a  spherically  symmetric 
electron  density  distribution,  and  that  the  effect  of  absorption  is  negligible.  As 
indicated  by  diffusion  theory,  the  electron  density  is  considered  to  have  a  Gaussian 
distribution  in  the  radial  direction. 

Initial  calculations  by  the  Haselgrove  procedure  were  presented  at  the  Four¬ 
teenth  Symposium  of  the  Electromagnetic  Wave  Propagation  Committee  of  AGARD 
in  Norway  (Rosenberg  et  al,  1968).  At  that  time  difficulties  were  encountered  with 
the  determination  of  the  path  for  the  ordinary  ray  at  small  input  parameters  where 
the  ray  reaches  the  critical  surface  and  a  cusp  or  "spitze"  phenomenon  occurs 
(Poeverlein,  1950).  The  scattering  amplitude  for  the  ordinary  ray  exhibited  very 
low  values  in  the  backscatter  region.  Inasmuch  as  these  results  did  not  appear 
trustworthy,  they  were  left  incomplete.  However,  since  that  time  the  computer 
program  has  been  improved  to  the  effect  that  the  ordinary  ray  can  be  calculated 
accurately  in  the  region  of  small  impact  parameters;  the  cross  section  no  longer 
exhibits  very  low  values  in  the  backscatter  region. 


2.  ANALYSIS 

A  ray  initially  a  distance  b  from  the  x  axis,  or  axis  of  propagation,  is  deflected 
through  an  angle  8  as  it  emerges  from  the  cloud  (Figure  1).  The  Haselgrove  equa¬ 
tions  for  the  ray  paths  are 


-Tg-(fi  cos  a  +  -g^sin  a). 


=  -tj(^»  sin  a  - cos  a ).  (2) 

i  (3) 

Here  /i  is  the  index  of  refraction,  c  the  velocity  of  light,  and  a  the  angle  the  wave 
normal  makes  with  the  direction  nf  propagation.  Since  the  constant  factor  c  has 
no  effect  upon  the  geometry  of  the  ray,  it  may  be  absorbed  conveniently  in  the  time 
factor  t.  The  index  of  refraction  is  given  by  the  Appleton-Hartree  equation,  which 
for  the  case  of  no  absorption  has  the  form  (Kelso,  1964a) 


tj*'*'*; 


2X(1-X) 


2(1 -X)  -  Yt2  ±  [Yt4  +  4Yl2  (1-X)2]1/2  * 


where  X  = - »,  «  and  co  are  the  plasma  and  incident  frequence  <a,  Y™,  and  YT  are 

WZ  n  1 

the  transverse  and  longitudinal  components  of  the  gyromagnetic  ratio  Y,  ari^  6  is 

the  angle  between  the  direction  of  propagation  and  the  earth's  magnetic,  field 

The  (+)  sign  refers  to  the  ordinary  ray  and  the  (-)  sign  to  the  extraordi  iary  ray 

Because  of  the  form  of  Eq.  (4),  the  scattering  phenomena  will  have  a  pt  vied  of  180° 

with  results  at  -8  the  same  as  those  at  180°  -  0  . 

P  P 

If  the  incident  rays  are  confined  to  the  magnetic  meridian  plane,  that  .  t  ie 

plane  containing  the  direction  of  propagation  and  the  magnetic  field  direction,  then 
no  lateral  deviations  occur  (Budden,  1966).  However,  for  rays  outside  of  thic-  plane, 
lateral  deviations  do  occur,  but  numerical  solutions  indicate  these  to  be  small 
(Kelso,  1964b).  The  effect  of  lateral  deviations  will  not  be  considered  here. 

For  a  Gaussian  distribution,  the  electron  density  n  has  the  form 


n  =  nQ  exp(-r  /&)  , 


where  nQ  is  the  center  point  density  and  a  the  Gaussian  radius.  If  we  define  a 
parameter  rc  by 


rclrcl  =  to(i*)  ’ 


where  nc  is  the  critical  plasma  density  for  the  ordinary  ray  (X  =  1 ),  the  electron 
density  may  then  be  written  as 


IT  =  exp  t(rc|rcl  "  r2,/a2l  ’ 


a  plot  of  which  is  shown  in  Figure  2.  We  note  that  for  overdense  plasmas,  nc 

occurs  at  a  radius  equal  to  r  ,  that  is,  r  is  the  critical  radius  at  which  the  local 

c  c 

density  is  critical.  For  plasmas  underdense  at  the  center,  rc  is  negative  but  sig¬ 
nificant  refraction  still  may  occur.  For  convenience,  we  shall  use  the  descriptive 
term  critical  radius  for  both  positive  and  negative  r£ .  Referring  to  Figure  3, 
which  gives  ray  diagrams  for  both  overdense  and  underdense  spheres,  we  note  that 
a  series  of  rays  coming  from  the  left  are  refracted  by  different  amounts  depending 
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on  the  impact  parameter  b,  the  initial  z  distance  above  the  axis  of  symmetry. 

For  large  values  of  b  there  is  little  or  no  refraction.  For  a  ray  at  small  b, 
bending  occurs  to  a  final  angle  ranging  from  0°  through  90°  to  180°  as  the  impact 
parameter  decreases  to  zero,  or  from  the  lower  side  from  0°  to  -90°  to  -180°,  as 
the  impact  parameter  increases  to  zero.  Note  that  the  direction  -180°  is  identical 
to  +180°,  both  indicating  complete  backscatter. 

In  an  underdense  plasma,  the  refraction  increases  from  0°  to  some  maximum 
angle,  say  90°,  returning  to  0  as  the  impact  parameter  reaches  zero.  Thus,  com¬ 
plete  backscatter  never  occurs  and  therefore  all  energy  is  scattered  in  the  forward 
direction. 

From  the  solution  of  the  Haselgrove  equations,  which  yield  the  increment  in 
scattering  angle  A0  corresponding  to  an  increment  Ab  in  the  impact  parameter, 
the  cross  section  a (6)  is  calculated  from 


<r(0) 


4lt  bAb 
sin  0  Ad 


(8) 


The  scattering  angle  and  the  cross  section  may  then  be  presented  as  a  function  of 
critical  radius,  frequency,  magnetic  field  direction,  and  polarization  (ordinary  or 
extraordinary  ray). 


3.  METHOD  OF  CALCULATION 

The  computer  solution  of  the  ray  tracing  problem  is  based  on  the  Haselgrove 
set  of  differential  equations,  which  fixes  ray  and  phase  angle  from  local  density  and 
density  derivatives.  With  the  assumed  Gaussian  distribution,  ordinary  and  extra¬ 
ordinary  rays  were  traced  for  a  matrix  of  two  frequencies,  four  magnetic  field 
angles,  and  seven  values  of  critical  radius,  with  a  fixed  Gaussian  radius  of  1  kilo¬ 
meter.  The  earth's  magnetic  field  Hg  was  taken  as  0.36  gauss. 

Rays  were  launched  in  the  x  direction  from  a  surface  8  km  from  the  center  of 
the  plasma,  at  a  series  of  21  z  distances  (impact  parameters)  above  and  below  the 
center  (Figure  3),  ranging  from  -4  to  +4  km;  each  ray  trace  ended  when  the  ray 
passed  out  through  the  8-km  surface,  noting  the  final  ray  angle  with  respect  to  the 
incident  ray  direction.  Steps  in  each  ray  trace  were  limited  to  a  maximum  time 
increment  chosen  to  limit  the  change  in  refractive  index  to  less  than  20  percent  and 
further,  to  limit  the  change  in  ray  angle  to  less  than  a  fixed  increment  d0  (1°,  0.5°, 
0.25°,  and  so  on). 


At  each  step,  a  Eunge-Kutta  approximation  to  the  Haselgrove  equations  was 
made  to  determine  increments  in  x,  z,  and  phase  angle,  which  were  then  used  to 
determine  the  effective  refractive  index  at  the  new  x,z  position. 

A  given  ray  tracing  was  repeated  with  successively  smaller  fixed  increment 
d@  values  until  successive  estimates  of  the  final  ray  angle  agreed  to  within  0.1°. 
Typically  100  to  600  steps  were  required,  and  2°  and  1°  fixed  increment  values 
agreed  to  0.1°. 

Cusping  of  the  ordinary  ray  was  indicated  by  the  phase  angle  approaching 
the  magnetic  field  direction  and  the  ray  angle  becoming  perpendicular  to  it.  When 
this  condition  was  encountered  to  within  0.3d6,  the  phase  angle  was  stepped 
"across"  the  magnetic  field,  reversing  the  ray  angle  automatically  by  180°,  the  ray 
tracing  being  continued  in  the  customary  manner.  The  region  in  which  cusping 
occurred  was  limited,  that  is,  if  the  ray  were  at  some  distance  from  the  critical 
surface  when  it  reached  the  magnetic  field  angle,  cusping  did  not  occur. 

The  axis  for  which  an  incident  ray  is  scattered  through  an  angle  of  180°  (over- 
dense  sphere)  or  0°  (underdense  sphere)  does  not  coincide  with  the  axis  of  sym¬ 
metry  of  the  cloud  because  of  the  asymmetry  introduced  by  the  magnetic  field.  The 
displacement  of  the  scattered  ray  was  obtained  from  interpolation  of  the  impact 
parameter  values  giving  both  positive  and  negative  scattering. 

Thus  the  output  of  the  ray  tracing  program  provides  scatter  angle  versus  im¬ 
pact  parameter  for  ordinary  and  extraordinary  rays  for  each  field  angle,  fre¬ 
quency,  and  critical  radius. 

4.  RESULTS  AND  DISCUSSION 

In  Figure  4,  some  typical  ray  paths  are  shown  for  zero  field  and  for  ordinary 
and  extraordinary  rays  at  a  frequency  of  3  MHz  and  a  field  angle  of  -45°.  The 
ordinary  and  extraordinary  rays  show  opposite  asymmetries  at  a  given  field  angle; 
for  example,  at  ®p  =  -45°  the  ordinary  ray  is  displaced  upward,  whereas  the  extra¬ 
ordinary  ray  is  displaced  downward.  At  a  field  angle  of  45°  (not  shown),  the  dis¬ 
placements  are  reversed,  indicating  a  cyclic  behavior.  For  an  underdense  sphere, 
the  effect  of  a  magnetic  field  is  shown  for  the  ordinary  ray  with  rc  =  -0. 50  and 
0p  =  -45°.  The  upward  displacement  of  the  rays  is  quite  significant,  leading  to  a 
large  negative  displacement  for  the  ray  scattered  in  the  forward  direction. 

The  variation  of  scattering  cross  section  with  scattering  angle  for  zero  field 
with  critical  radius  as  parameter  is  presented  in  Figure  5.  The  corresponding  re¬ 
sults  for  the  magnetic  field  case  are  presented  in  Figures  6  through  9  for  frequen¬ 
cies  of  3  and  30  MHz  and  for  field  angles  of  0°,  ±45°,  and  90°.  The  cross  section 
has  been  expressed  in  terms  of  effective  radius  re  corresponding  to  the  radius  of  a 
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perfectly  conducting  sphere  with  the  same  cross  section  as  the  given  sphere. 

Since  there  is  little  ray  penetration  into  a  highly  overdense  plasma  (large  rc),  and 
conversely,  deep  penetration  into  a  slightly  overdense  or  underdense  plasma  (rc 
small  or  negative),  it  will  be  convenient  to  refer  to  these  as  "hard"  and  "soft" 
spheres,  respectively. 

For  zero  field,  the  cross  section  for  hard  spheres  is  high  in  the  backscatter 
region,  with  very  little  change  over  a  wide  range  of  scattering  angles.  As  the  crit¬ 
ical  radius  decreases,  the  backscatter  cross  section  decreases  rapidly  along  with 
a  decrease  in  the  range  of  substantially  constant  cross  section.  As  we  approach 
the  critical  density,  rc  =  1,  the  cross  section  becomes  negligible  in  the  region 
from  90°  to  180°.  This  result  is  indicated  by  a  vertical  line  at  B  =  90°.  Forunder- 
dense  spheres,  the  range  of  scattering  angles  is  confined  to  the  forward  region,  the 
cross  section  becoming  large  in  the  neighborhood  of  the  maximum  scattering  angle. 
This  effect,  attributable  to  a  small  change  in  scattering  angle  with  impact  param¬ 
eter  in  the  region,  has  been  indicated  by  a  vertical  line  at  maximum  scattering.  We 
note  that  the  curves  for  different  critical  radii  converge  in  the  region  of  forward 
scatter,  coalescence  occurring  at  smaller  angles  as  the  critical  radius  increases. 
This  behavior  is  derivable  from  the  nature  of  forward  scattering  for  which  the  im¬ 
pact  parameter  is  large,  with  ray  deflection  occurring  in  the  outer  region  where 
t’»e  electron  density,  together  with  its  gradient,  is  small.  Scattering,  therefore, 
depends  principally  upon  the  impact  parameter  with  little  dependence  upon  the  crit¬ 
ical  radius.  The  effect  of  critical  radius,  however,  is  not  negligible  when  rp  is 
large  (see  Figure  2),  resulting  in  the  delayed  convergence. 

In  general,  the  effect  of  a  magnetic  field  is  quite  small  at  a  frequency  of  30  MHz, 
becoming  appreciable  at  3  MHz.  The  ordinary  ray  exhibits  a  minimum  at  a  field 
angle  of  90°;  the  minimum  becomes  more  pronounced  as  the  frequency  is  reduced 
from  30  to  3  MHz.  A  wavelike  pattern  is  evident  in  several  of  the  curves,  particu¬ 
larly  for  the  softer  spheres  at  the  lower  frequency.  It  is  difficult  to  find  a  simple 
explanation  for  this  behavior  because  of  the  complexity  of  the  computing  process. 
However,  a  careful  examination  of  the  data  indicates  that  the  behavior  is  not 
spurious. 

The  curves  for  the  extraordinary  ray  are  quite  smooth,  bearing  a  closer  re¬ 
semblance  to  the  zero  field  results  than  do  those  for  the  ordinary  ray.  At  a  fre¬ 
quency  of  30  MHz,  the  results  are  very  close  to  those  for  the  zero  field.  However, 
the  critically  dense  sphere  now  exhibits  definite  overdense  characteristics,  that  is, 
a  nonnegligible  cross  section  in  the  backscatter  region  and  a  maximum  scattering 
angle  of  180°.  As  the  frequency  is  decreased  to  3  MHz  the  backscatter  cross  sec¬ 
tion  rises  substantially  but  the  form  of  the  curves  remains  unaltered.  In  addition, 
the  range  of  overdense  characteristics  now  extends  to  an  underdense  sphere  with  a 
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critical  radius  =  -0.5.  It  thus  appears  that  the  principal  effect  of  the  field  with 
respect  to  the  extraordinary  ray  is  that  of  increasing  the  effective  critical  radius 
of  the  cloud,  giving  it  the  characteristics  of  a  hard  sphere. 

To  show  more  pointedly  the  relation  between  cross  section  and  scattering 
angle,  polar  diagrams  are  presented  in  Figures  10  through  12  for  several  field 
angles.  To  avoid  negative  values,  the  length  of  the  radius  vector  has  been  taken  as 
1  +  log  rfi.  Where  the  scattering  is  confined  to  the  forward  direction  the  maximum 
deflection  angle  has  been  indicated  by  a  straight  line  from  the  origin.  Because  of 
the  closeness  of  the  curves  in  the  forward  direction,  the  corresponding  values  of  rc 
have  been  marked  with  respect  to  the  maximum  deflection  lines. 

As  indicated  previously,  the  calculation  of  the  cross  section  requires  the  dis¬ 
placement  6  from  the  axis  of  symmetry  of  the  cloud,  of  the  ray  scattered  through 
180°  or  0°.  To  show  the  characteristics  of  the  displacement,  we  have  plotted  it  as 
a  function  of  field  angle  for  several  values  of  critical  radius  in  Figure  13.  Both  the 
ordinary  and  extraordinary  rays  show  a  definite  cyclic  behavior  with  respect  to  0p. 
The  magnitude  of  the  displacement  is  a  maximum  near  45°  and  135°  and  a  minimum 
in  the  neighborhood  of  0°  and  90°.  For  the  ordinary  ray  0  is  negative  for  0p  in  the 
range  0°,  90°  and  positive  in  the  range  90°,  180°,  whereas  for  the  extraordinary 
ray  the  sign  is  reversed.  We  may  anticipate  this  behavior  since,  for  example  (Fig¬ 
ure  4)  the  ordinary  ray  at  0p  =  -45°  range  shows  a  downward  bias  near  the  axis  of 
symmetry  of  the  cloud,  requiring  a  positive  correction.  The  displacement  is  quite 
small  for  the  extraordinary  ray  at  both  frequencies,  but  it  shows  a  moderate  in¬ 
crease  for  the  ordinary  ray  when  the  frequency  is  decreased  from  30  to  3  MHz.  To 
show  the  behavior  of  the  displacement  over  the  entire  range  of  critical  radius,  we 
have  plotted  0  as  a  function  of  rc  for  a  field  angle  of  45°  (Figure  14).  The  dis¬ 
placement  is  a  maximum  at  rc  =  0  for  the  ordinary  ray,  whereas  for  the  extra¬ 
ordinary  ray  a  maximum  occurs  at  rc  =  -0. 5  at  3  MHz.  For  the  extraordinary  ray 
at  30  MHz  and  similarly  for  the  other  field  angles  (not  shown),  the  displacements 
are  too  smaU  for  any  great  significance  in  the  patterns. 

To  exhibit  explicitly  the  effect  of  field  angle  upon  the  cross  section,  we  have 
plotted  re  against  0p  with  rc  as  parameter  (Figure  15).  For  reference  the  cor¬ 
responding  values  for  zero  field  are  indicated  in  Figure  15a  (ordinary  ray).  The 
ordinary  ray  shows  a  definite  maximum  at  0p  =  90°  for  both  frequencies  at  back- 
scatter.  Further,  to  indicate  the  behavior  at  another  scattering  angle  we  have  in¬ 
cluded  the  results  for  90°  scatter  at  3  MHz.  The  pattern  is  quite  different,  with  a 
minimum  near  0p  =  45°  and  a  maximum  near  135°.  The  extraordinary  ray  shows 
a  similar  pattern  except  that  a  minimum  is  replaced  by  a  maximum  and  vice  versa, 
with  variation  not  as  great.  In  fact,  at  30  MHz  the  rc  =  0  curve  shows  no  depend¬ 
ence  upon  field  angle. 
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The  dependence  of  the  maximum  scattering  angle  for  underdense  spheres  upon 
field  angle  and  sphere  hardness  is  shown  in  Figure  16  where  ®iaax  is  plotted  as  a 
function  of  8p  for  several  values  of  r£.  The  ordinary  ray  generally  shows  a  weak 
dependence  upon  field  angle  except  at  3  MHz  and  the  higher  values  of  rc  where  the 
dependence  can  be  considered  as  weak  to  moderate.  For  the  extraordinary  ray, 
there  is  virtually  no  dependence  at  30  MHz  and  a  rather  weak  variation  at  3  MHz. 

We  have  noted  heretofore  that  the  backscatter  cross  section  for  the  extraordi¬ 
nary  ray  exhibits  a  strong  similarity  to  that  obtained  for  a  hard  sphere  when  no 
field  is  present,  that  is,  high  values  of  the  cross  section  and  a  large  range  of  back¬ 
scatter  angles  over  which  the  cross  section  is  substantially  constant.  Furthermore, 
the  dependence  of  the  cross  section  upon  field  direction  is  rather  weak.  It  is, 
therefore,  natural  to  inquire  whether  the  backscatter  data  for  the  extraordinary  ray 
can  be  represented  in  terms  of  the  critical  radius  for  no  field  and  the  gyromagnetic 
ratio.  Because  of  the  small  effect  at  a  frequency  of  30  MHz,  we  confine  our  anal¬ 
ysis  to  3  MHz,  where  Y  =  1/3.  Since  the  critical  point  for  the  extraordinary  ray 
occurs  at  X  =  1  -  Y,  a  possible  relation  for  a  Gaussian  distribution  is 


r  2  =  r  2  -  ln(l-Y) 
ce  c 


=  r  ^  +  0.40  for  Y  =  1/3  , 


(9) 


where  r  is  the  effective  critical  radius  which  yields  the  same  cross  section  at 

2  2 

zero  field  as  obtained  with  the  given  field.  A  plot  of  r  against  r  (Figure  17) 

2  2  c 

shows  that  rce  is  indeed  a  linear  function  of  rc  except  at  small  values  of  rce 
where  we  do  not  expect  the  hard  sphere  concept  to  be  accurate.  The  dependence  of 
rce  upon  field  angle  shows  up  in  the  deviation  in  both  intercept  and  slope  of  the 
straight  line  representing  the  data  from  the  45°  line  given  by  Eq.  (9).  The  results 
for  45°  and  -45°  are  virtually  the  same  and  fairly  close  to  the  45°  line,  with  some¬ 
what  larger  displacements  for  the  0°  and  90°  angles.  Attempts  to  represent  the 
angular  dependence  in  terms  of  an  effective  average  angle  for  the  ray  path  have  not 
proved  fruitful.  It  is  possible,  however,  to  represent  the  data  in  a  simple  manner 
in  terms  of  a  modified  form  of  Eq.  (9),  that  is, 


rce2  =  AK2  "  ln(1-Y)J  +  B> 


(10) 


where  A  and  B  are  functions  of  8  .  Evaluation  of  the  A  and  B  terms  and  using 

P 

a  parabola  in  sin  8  to  express  field  angle  dependence  yields 


A  =  1.055  -  0.090  sinV 

P 

(ID 

B  =  0.15  -  0.23  sin2fi  . 

P 

Although  calculations  at  higher  values  of  Y  are  not  available,  we  may  expect 
Eq.  (10)  to  be  a  valid  representation  of  the  data.  In  fact,  in  view  of  the  augmented 
hardness  of  the  sphere  at  larger  Y,  the  A  and  B  coefficients  should  exhibit  a 
weaker  dependence  upon  field  angle  as  Y  increases.  Since  the  effect  of  field  angle 
is  not  large,  it  is  .’seful  to  note  the  degree  of  error  incurred  in  the  cross  section 
by  use  of  the  simpLfied  Eq.  (9)  in  place  of  the  data.  Although  the  deviation  in  rcg 
is  greatest  at  large  values  of  r  ,  the  accuracy  is  best  here  because  of  the  weak 
dependence  of  cross  section  upon  critical  radius  when  it  is  large.  For  example  at 
rc  =  1.41  where  log  r0  =  0.14  at  6p  =  0°,  0.12  at  45°  and  0.08  at  90°,  Eq.  (9)  yields 
0.11;  at  rc  =  0,  where  log  rg  =  -  0.35,  at  0°,  -0.42  at  45°  and  -0.51  at  90°,  Eq.  (9) 
yields  -0.47.  The  cross  section  obtained  from  Eq.  (9)  may  therefore  be  used  as  a 
reasonable  average  estimate  when  field  direction  effects  are  ignored. 
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Figure  3.  Ray  Diagrams  for  Different  Values  of  Critical 
Radius,  30  MHz  Extraordinary  Ray,  0p  =  0°. 

(a)  r£  =  1.0;  (b)  rc  =  0;  (c)  rc  =  -1.0;  (d)  rc  =  -1.41 


Figure  4.  Typical  Ray  Paths  Through  Gas  Cloud. 

(a)  r£  =  1.0,  no  field; 

(b)  rc  =  0.5,  3  MHz,  extraordinary  ray,  8p  =  -45°; 

(c)  rc  =  1.0,  3  MHz,  ordinary  ray,  8p  =  -45°; 

(d)  r£  =  -0.5,  3  MHz,  ordinary  ray,  8p  =  -45° 
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ordinary  ray,  30  MHz;  (b)  ordinary  ray,  3  MHz;  (c)  extraordinary  ray,  30  MHz;  (d)  extraordinary  ray,  3  MHz 
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Figure  10.  Polar  Diagram  of  1  +  log  r£  as  a  Function  of 
Scattering  Angle  for  8p  =  0. 

(a)  ordinary  ray,  30  MHz;  (b)  extraordinary  ray,  30  MHz 


(b) 


(d) 


Figure  11.  Polar  Diagram  of  1  +  log  rg  as  a  Function  of 
Scattering  Angle  for  8p  =  45°. 

(a)  ordinary  ray,  30  MHz  (b)  ordinary  ray,  3  MHz 

(c)  extraordinary  ray,  30  MHz  (d)  extraordinary  ray,  3  MHz 


(b) 


.  ,  +  log  r  as  a  Function  of 

Figure  12.  Polar  ^agr^^ 

Scattering  Angle  10  up 

_  in  MHz 


Scattering  Angle  for  8  -  •  ^  ordinary  ray,  3  MHz 

(a)  ordinary  ray,  30  MHJ  MHz  (d)  extraordinary  ray, 
(c)  extraordinary  ray,  » 
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Figure  14.  Displacement  as  a  Function  of  Critical  Radius 
With  Frequency  as  Parameter;  0p  -  45 
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